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SUMMARY 


The potential of the linearized flow for a lifting-line propeller 
of arbitrary circulation distribution at subsonic advance is derived in 
cylindrical harmonics . From the potential the induced angle at the 
lifting line is obtained. The series expression for the induced angle 
is divergent as is to be expected for lifting lines in supersonic flow, 
but this divergence is removed when the supersonic lifting-line induc- 
tion is removed. A phenomenon resembling resonance in vibrating systems 
introduces itself in that one term in the series becomes very large com- 
pared to the others. The main consequence of this phenomenon is that 
the lift distribution cannot be arbitrarily prescribed; on the other 
hand, the inverse problem, in which the propeller geometry is given, 
is acceptable. 


INTRODUCTION 


In order to calculate the induced angle of attack of the propeller 
blade sections for compressible flow, the velocity potential for the 
complete flow field must be derived. Inasmuch as the equation of motion 
changes from an elliptic type to a hyperbolic type at the radius where 
the propeller blade, or the blade prolonged, is at Mach number one, the 
disturbance potential turns out to be of mixed elliptic and hyperbolic 
character. With incompressible flow, the distant boundary condition 
is certainly that the disturbance potential must be zero far out radially 
and far upstream. However, with compressible flow, the possibility of 
propagation of disturbances far ahead, especially in a closed wind tun- 
nel, must be admitted. Therefore, the one problem becomes two; the 
elliptic field, which is essentially incompressible, and the hyperbolic 
are determined by different physical conditions. 

The velocity potential is obtained by superposing suitable solutions 
of the compressible equation of motion on the known far -wake velocity 
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potential, which is also a solution of the compressible equation of 
motion. This far-wake potential extends only downstream of the pro- 
peller and provides the jump in potential at the trailing vortex sur- 
faces and the downstream feature that the far wake must be the same as 
for incompressible flow except for random wave motion or noise. There- 
fore, the problem reduces to superposing the elliptic and hyperbolic 
fields in such a manner that certain physical conditions are met at 
the propeller plane where the far-wake potential is cut off, an obvious 
condition being continuity of the velocity vector. The cut-off of the 
far-wake potential at the propeller plane creates the lifting lines 
there, since this termination of the surfaces of a potential discon- 
tinuity is equivalent to a lifting line. 

In the differentiation of the potential at the lifting line, a 
difficulty is to be anticipated in that a lifting line has infinite 
wave drag in supersonic flow, and this condition is part of the theory 
inasmuch as no restrictions are placed on tip Mach number. Some of 
the hyperbolic induction attributed to the lifting line must therefore 
be separated. This unwelcome induction is found in the hyperbolic 
field but is absent in the elliptic field, as would be expected. 

For the most part, only the case of propeller operation in a closed 
circular wind tunnel is considered, although the way the potential may 
be obtained in free air is indicated. The only reason for emphasizing 
the tunnel is that, where series occur in the tunnel theory, integrals 
occur for free air. The series seem more amenable to investigation 
than the integrals, and it is believed that free air is the limiting 
case in which the tunnel diameter approaches infinity. 


SYMBOLS 


B 


number of blades 


V 


M 


advance Mach number, v/a 
advance velocity 


R.P. 


real part 


a 


velocity of sound 



t 


time 


axial disturbance velocity, positive downstream 
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u t tangential disturbance velocity, positive in direction of 

propeller rotation 

u z> u D u a disturbance velocities in positive directions of coordi- 
nates z, r, and a, respectively 

y nBj argument of jth extreme value of Bessel function of order nB 


x,p,a 


z,r,a 

r(p) 

0 


cp,f 

Oi 


dimensionless cylindrical coordinates (eq. (4)) 

cylindrical coordinates 

circulation at p 

disturbance velocity potential 

complementary parts of velocity potential 0 

induced angle of attack at blade 


3 = \j± - M 2 


£ helical angular coordinate (eq. (7)) 

p a.’ p T’ p S particular values of p for propeller tip, tunnel wall, 
and sonic radius, respectively 

o> angular velocity, radians per unit time 

Subscripts : 


i incompressible 

c compressible 


Partial differentiation is indicated by subscripts . 


VELOCITY POTENTIAL FOR FLOW IN A CLOSED CIRCULAR WIND TUNNEL 

For a linearized theory the disturbance velocity potential 0 
must satisfy the equation 

a 2 ” ^ r 2 r 2 + r£ ^r 2 + ^2 + ^z 2 z 2 


( 1 ) 
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in which Z2,r2,a 2 are cylindrical coordinates fixed in the distant 
fluid, t is the time, and a the velocity of sound. In coordinates 
which rotate and advance with the propeller, the flow is steady after 
the propeller has been in operation for a sufficiently long time under 
constant conditions . Let z-^,r-^,aq be the coordinates in the moving 
system with zq along the propeller axis of rotation and with positive 
directions downstream, outward from the axis, and counter to the pro- 
peller rotation, respectively (see fig. l). If V is the velocity of 
advance of the propeller and cd the angular velocity, substitution of 
the transformations 


zi = z 2 + Vt 


a l = a 2 + ^ 


r l ~ r 2 


( 2 ) 


into equation (l) leads to the steady-flow equation 


Vf 

ZtZi \- „2 J ' rr i r i ' rq ''r-j 




— 0 (1 

2 


ou?r n 2 ' 


- 5 ■ 0 


(3) 


This equation has a more compact form in a third coordinate system x,p,a 
differing from the Zq,rq,a.q system only in the length scale; that is, 


o>zq cor-i 

x = p = — 1 a = an 

V V 


(M 


In this new coordinate system equation (3) becomes 


- M 2 ) + 0 pp + i 0 p + ~ 0^(1 - P 2 M 2 ) - 20 xa M 2 = 0 (5) 


In the helical coordinates of Goldstein and Reissner (refs. 1 and 2, 
respectively) equation (5) becomes 


0xxP 2 + 0pp + ^ 0p + 0^(d + P~ 2 ) _ 20 £x = 0 


( 6 ) 


MCA TN 2983 


5 


through the transformation 


£ = a - x (7) 

In the far wake, the oblique system used in equation (6) has advantages 
because of the helical nature of the flow. However, except for far-wake 
considerations, equation ( 5 ) is preferred because of its orthogonal 
coordinate system. In equation (6) if M = 1, then (3 = 0 and the 
equation suffers the degeneration typical of the linearized equations 
of flow in that the coefficient of the (f ^ term is no longer repre- 
sented by p with sufficient accuracy. Therefore the advance Mach 
number must not be too close to one . 

Equations ( 5 ) and (6) may be shown to change from an elliptic type 
to a hyperbolic type on a cylinder concentric with the propeller axis 
of rotation with a radius such that the corresponding p is 


(8) 


where 0 = vl - M^. As one would expect, equation (8) defines the 
cylinder on which the resultant velocity of the undisturbed stream in 
the moving coordinates is equal to the speed of sound. Because of the 
change in type of the equation of motion, mixed elliptic and hyperbolic 
components will generally occur in the complete solution. 

The potential given by Reissner in reference 2 for the far wake 
is already in a usable form for compressible flow because it satisfies 
equation ( 5 ) or (6). This potential is for a propeller of arbitrary 
circulation distribution. Adapted to the case where a circular wind 
tunnel is present, Reissner' s potential 0^- for the far wake takes 
the form 





BT(p) 

2k 


z + 


00 



-2(-l) n 


nB 


h n (p)sin nB£ 


(9) 
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where 


h n (p) = -K(p) / I'(|)6~d / I ' ( | ) | 7— d| - 


'0 


d| 


I(p)k'(p t ) r p a d£ 

'0 d £ 


I'(P T ) 


I(p) / a K'(|)| ^ d| 

Jp 


(p < p a ) ( 10a) 


h n(p) 


-K(p)l'(p T ) + l(p)K'(p T ) 
I'(P T ) 



rdhlai 


(p a < P < Pt) 

(10b) 


In these equations p a and p>j are the values of p for the 

propeller tip and the tunnel radius, respectively, B is the number of 
blades, f = T(p) the circulation distribution depending only on p, 
and K(p) and l(p) are abbreviations for K nB (nBp) and I nB (nBp), 

respectively, which are Bessel functions of imaginary argument. The 
primes on K’(p) and I'(p) denote differentiation with respect to p, 
not with respect to the argument nBp. The angle £ is zero halfway 

between vortex sheets from adjacent blades and is limited to - £ < £ < §■ 

"Vrr B -D 

with the vortex sheets at ± S-, ± — , .... 

•D R 


It is immediately apparent that the Reissner potential, independent 
of x, is a solution of equation (6) and hence of equation ( 5 )? because 
for a potential independent of x, equation (5) is Laplace's equation 
in the p,£ system. Therefore, the Reissner potential equation (9) is 
considered as the first component of the velocity potential for the pro- 
peller with compressible flow, but its region of application is only in 
back of, or downstream of, the propeller plane. Thus, the jump in (J) 
at the helical surfaces trailing from the blades is provided. Therefore, 
the flow divides naturally into two regions , one in front and the other 
in back of the propeller plane, both inside the tunnel wall. The rest 
of the problem consists in adding suitable solutions of equation (5) in 
front and in bad: of the propeller in order to take account of the 
physical conditions . 
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The appropriate solutions of equation ( 5 ) are 



as may be verified by direct substitution. The constants n and l are 
free. The tunnel-wall boundary condition requires that 


PTj/^nj 2 + (nB)2 ^ ' = < 12 ) 

where y^gj is the argument of the Bessel function in equation (ll) at 
the jth extreme value. From equations ( 8 ) and (l2) 


nj 



(13) 


where the plus sign is chosen to define l n j . 


Equation (ll) now becomes 


M-F “ 


B 



±Z„ -j+i 


nB 


nj Ps 2 y 


< 


>e iriEa< 


> 



s(-ilnk)+25_ 


S J 


(14) 


The notation p F refers to p in front or back of the propeller plane, 
B 

the F going with the upper sign in the braces, and the B with the 
lower. The meaning of the two functions, one above the other, in braces 
is that the upper function is used when is real, the lower, when 
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imaginary. It is convenient to use l n ^ to indicate the imaginary 
values. The upper function is called, for convenience, the elliptic 
solution because it dies out exponentially with x, and the lower is 
called the hyperbolic because it oscillates without dying out. This 
complication, two possibilities for function character, must be accepted 
because of the mixed flow. 

Two more physical conditions, altogether different, have been 
imposed in the choice of ± instead of + in equation (lA) for the two 
functions p in the braces . The origin of x has been set at the 
propeller plane with x positive downstream; for the elliptic functions 
the choice of signs is obviously that which makes the solution vani sh at 
x = ±00. The choice in the hyperbolic case is such that, since -li nk 

is positive, the solutions oscillate for given values of n and k 
with higher frequencies in front than in back of the propeller, in 
accordance with the physical observation in reference 3 that waves are 
crossed with higher frequency in front than in back, proceeding in the 
axial direction. 

The rest of the problem is to make superpositions of equation (l4) 
on the Reissner potential, which occupies only the region in back of 
the propeller, in such a way that the flow is continuous in the propeller 
plane off the lifting lines. 

Because of the discontinuous nature of the first term in equa- 
tion (9), the compressible fields which must be superposed on the two 
terms composing equation (9) must be determined separately. For this 
purpose, let 


= % + ^ + % + *0B + ^ (0 < x) (15) 


0 F = % + *0F + ^F (x<0) (16) 

in which i]/ R is the first, or discontinuous, function in Reissner 's 

potential (eq. (9 )), is the second, or continuous one, fp and cp^, 

B B 

are the solutions of the form of equation (l4), and \|r QF are degenerate 

B 

solutions of the form of equation (l4) with n = 0. Now the origin for 
the coordinate a must be decided; this is taken to be halfway between 
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A 


blades. The angles a and £ are now the same at x = 0, with the 
nearest blades (lifting lines) at a = £ = ± x = 0. The angles 

£ = ± ^ in the x,p,£ system are at the nearest trailing vortex sheets 
on each side of £ = 0. 

The compressible potentials ^ and ^ which accompany the dis- 
continuous Reissner function \|ig are now determined to provide con- 
tinuity at the propeller plane. From equation (9) 


% 


- _ X ) 


2k 


(17a) 


*R =R - P - IfX 

n=l j=l 


- g( ; l)tl r p)e lnB(a - x) (-i 


2 k 'rzv nB 1 nj°nB[ p p 


(-i) ( 17 b) 


where the second form is obtained from an expansion of a - x in a 
Fourier sine series and another expansion of r(p) in the form 


00 


F(p) - y Tnj'-'nB 
j=l 



(18) 

1 


in which r n j is real. This Fourier-Bessel type of expansion is used 

repeatedly herein. In particular, the arguments of the Bessel functions 
are always adjusted through the frequency coefficient y^Bj/PT "to make 

the extremes fall at the tunnel wall prp ; in which case r n j is given 
by (see ref. p. I 7 M 
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The function which is limited to the region downstream of the 

propeller, causes a discontinuity in the axial velocity and in the 
potential at the propeller plane x = 0. In order to avoid differentia- 
tion of equation (17b), the discontinuity in axial velocity brought by 
can be obtained from equation (17a) directly and eliminated with a 

degenerate term from equation (ik) obtained by setting n = 0; the 
result is 



tOF " 


B 

in which Aqj 

is to be 

introduced by 

\|/Qp , and 


B 

by 



^ p ) e 


P T 




(20) 



BT(p)/ \ 

- x) 

+ 

^OB 



_dx 

2jt J 

x=0 

_8x_ 

x=0 

_8x_ 


xO 


( 21 ) 


From equations (17a) and. (20), equation (2l) becomes 


2* 2rt 4 = T ROm ^ 0 \P T / 2lt 0P T ^ °\ Pm ) 


j=l PP T 

By using equation (l8), the coefficients A qj become 


T 


( 22 ) 


1 PPt 

= - 5 — r 0j 


y 0j 


(23) 


The discontinuity in axial velocity brought by the Reissner term ^ 
is therefore canceled by superposing the potentials 
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I 



(24) 


Only the discontinuity in potential at x = 0 caused by 
remains. This discontinuity is canceled by requiring that 





(25) 


and, in order to prevent introduction of a discontinuity in axial 
velocity at x = 0, by requiring that 


ft ‘ ft % ] 


(26) 


x=0 


in which \|r R comes from equation ( 17 b), and \|r F is written in the form, 

B 

adopted from equation ( lk ) , 


r 


j=°° 

*_ = R.P. ^-Z^- ~ g( - 1>n 

B 2 ” Si 

k=l 


A njF 

B 


A nkF 

B 


'nB 


^nBj 




re 


T (y nBk ^ 

H *- p j 


inBa^ 


±Zn j +i ^l2 y ,x 


i (+ iZ nk + ^2' x 


( 27 ) 


in which A n j and A ^ “ay be complex; that is, 


A nj a nj + -*- a nj 


A nk a nk + ^ a nk 
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The range k = 1, 2, 3, . . . kQ denotes the hyperbolic terms, and 
j = j 0 > jo+1’ ^0+2’ • • • 00 denotes the elliptic. Here, denotes 

the last hyperbolic term in the summation and . Substituting 

equations (l 7 h) and ( 27 ) into equation ( 25 ) gives 




r f 




anjB + ia njB 

1 

>- -j 

anj F + ia njF 

► 

* 

^nk( ) 


a nl-B + la nkB_J 

I | 

a-nkF + ia nkF 
- 


and, similarly, substituting equation ( 27 ) into equation ( 2 6 ) gives 


-Z n j + i 25- 
P S 2 


r . . 

a njB + ia njB 

< 



iAink + ] 


a nkB + ia nkB 

V P SV 





+ ia njF 


+ ia nkF 


(29) 


Equations ( 28 ) and (29) each break into two equations corresponding 
to the real and imaginary parts which yield the following solutions for 
the elliptic field: 


a njF " a njB 


1 nB 

2 


z njPs 


2 r nj 


> 


a njF “ a njB ~ 2 r *hj 


-y 


and for the hyperbolic field 


a nkB ® 


a nkF 


_ 1 
" 2 



nB 


-iZ n kPg 




a nkB 


. 1 
" 2 




(30) 


( 31 ) 


r 


The coefficients given by equations (30) and (31) are now to be used in equation (27). Then 
according to equations (15) and (l6), the \(r part of the potential given by collecting equa- 
tions (17a), (2b), and (27) is, for the region in back, 


+ *013 + 


Bf(p) 

2k 


(a 



i- r, 


0j 






+ 


0 =°° 

CO ^0 / ill 

R.P. IVV I 2 *- 11 
& fel 0=j o "B 2 

k=l 


h< 


nj 


nk 


-nB 


! noP s 2 


i/1 + 


+ i 


nB 


J nB 


f y nBj 


>< 


■iinkPg^ 


Jv e inBa^ 


/ 1 

fy n 


JnB • 


nBk 
\ P T 




. / . , nB 
1 ^ 1 ^nk + — 

P S 


and in front, 

'•'OF 


(32) 


yp.i 


k=t 


r.p. B.57 y -g (-i) n J 

9-n- ~ - 

k=l 


r nj 


-nB 


l njPs 2 


'nB 


2 * SL f^ 0 2 


nk 


< 


if -1 + 


nB 


“i^nkPg^y 


J nB 


y nBj 
v P T 

f T 

y nBk 


inBa< 


. nB , 


1,_1 W“5V 


04 


( 33 ) 
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Equations (32) and (33) comprise the final results for the \|r part 
of the potential 0. 

Nov consider the cp part. The procedure is nearly the same, only 
this time the discontinuity in axial velocity will not be dealt vith 
separately because the cpp part of Reissner’s potential is a smooth 
Fourier development in a - x which can be differentiated. The con- 
tinuity conditions at the propeller plane x = 0 become 


IhiU + WxK, ■ M 


1 x=0 


(34) 



+ 

Scpg 


8cpp 

_Sx _ 

x=0 

8x _ 

x=0 

8x_ 


^ J x=0 


(35) 


in which cp^ is the second term in equation (9)- This term may he 
written, in the x,p,cx system, as 


„ ^ B V” -2(-l) n , s inB(a-x) , . * 

<PR = R - P ‘ 2^^ Mp)e (-1) 

n=l nB 


(36) 


or, alternatively, 


„ „ B V" V" -2(-l) n v T ^nBj \ inB(a-x), ^ 


(37) 


where 


>n(p> * Zi VjWlT 1 

j=l \ m t 


(38) 


and h n j is real. 


I 


The 


cpp are obtained from equation (27) with * replaced by cp and some new coeffi- 
B 


cients substituted for A 


nj 


, fSt * Part ° f thG P° tential now becomes, with the B ni determined from equations ( 34 ) 
and (35) in a manner similar to that for A n j, J ^ w 


q>R + cpfi = R.P. ~ - 2 ^ l)n h n (p)e irlB(a - x) (-i) + 


R.P. - 2 - 


J=00 

k=kQ 


-2(-l) n l 


2,1 £1 J=r 0 112 2 

k=l 


nj 


4 ik 


-nB 


>< 


Z njP 2 


i/1 + 


+ i 


nB 


p) 


x 


‘ 1! njS 


\ 

j fy^Bk \ 

) 

p ) 


>e inBa J 


7 .nB , 


{ i!nk+ S 


> 


|X 


( 39 ) 


J=00 

k=i 



V 


^ = =s£=i£i< 

^ fee M 0 »B 2 s 

k=l 


-nB 


* 4 ik 




i/ -1 + 




— i 




nB 


-il nk 3 2 > 


r e 


inBa «! 


'nB 


f y nBk 

Pip 


ln j + W 


i ( “iink+T^j | x 




p S 


(40) 


Equations ( 39 ) and- (^-O) added to equations ( 32 ) and ( 33 ) determine the velocity potential 
for the propeller in accordance with equations (15) and (l6). 
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PRANDTL INDUCED VELOCITIES AT THE LIFTING LINE 


The Prandtl induced velocities, in the sense used herein, are analo- 
gous to those of the conventional lifting-line wing and propeller theo- 
ries of Prandtl and Goldstein. More specifically, the Prandtl induced 
velocities are those which must be presumed to be already present at 
the position of the blades in order to assume that the airfoils retain 
locally their same performance as in two-dimensional flow. With incom- 
pressible flow, it was verified experimentally that the induction could 
be determined with sufficient accuracy by assuming the wing or blade to 
shrink into a lifting line and then applying the Biot-Savart law over 
the trailing vortices excluding the lifting line. This simplification 
of the intrinsic general three-dimensional problem has long since 
proved its worth. Therefore it is natural to seek extensions in 
approaching the compressible problem. On the other hand, it is dan- 
gerous to rely on experience with incompressible flow for guidance in 
supersonic matters . Since the present problem concerns a mixed subsonic 
and supersonic flow, it will be expected that questions will arise which 
will have to be settled by considerations of both subsonic and super- 
sonic flow. Finally, for those who are more interested in the general 
three-dimensional problem rather than an engineering concept, it is 
remarked that the theory is still applicable to a lifting-surface 
problem when those steps are omitted which specialize to the lifting 
line . 

The induced velocities, or disturbance velocities, are given by 
differentiations of 0, but the algebraic sign is still open to a choice 
which is made so that a forward- thrusting propeller pushes the air rear- 
ward in the far wake. If the propeller is assumed to produce positive 
thrust when r(p) is positive, then the velocities in the positive 
directions of the coordinates must be 


which may be checked by means of the Reissner \|/p term, which produces 
the basic disturbance velocity in the far wake, as follows: 


Uz i = ~^ z i Ur i = ~^ r i X = ' ^ 


(4l) 



Since 


A 
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is seen to be negative, but since the propeller rotates in the 
negative exp direction, is a velocity in the same direction as 

the rotation, as it should be; on the other hand, u z -^ is positive, 


but since zq is positive downstream, 
also as it should be . 


u 9 


is positive downstream, 


In propeller theory, it is customary to consider an axial veloc- 
ity u a and a tangential velocity u^, which are positive downstream 
and in the direction of propeller rotation. Since u z -^ is positive 
downstream, but u^ is positive against the propeller rotation, 

u a = u z and u-t = -u . Therefore, from equations (Al) and (A) 

"1 ^1 


0) ri 

- - 0x 

(42) 

“ 1 gL 
v p ycL 

(43) 


determine the conventional axial and tangential velocities from the 
velocity potential in the x,p,a system. 

It is immaterial whether the differentiations in equations (42) 
and (45) are applied to or 0-g when the differentiation, which 

is not straightforward, is made at the- lifting lines. Suppose 0-g is 

used and consider first the \|r part of (f> (eq. (32))- For the \|r x 

differentiation, the second term in equation (32) cancels one-half of 

the first term. If it is permissible to differentiate the remaining 

term, then at x = 0 and a = which means- the lifting line, 

B 


dx 


(^R + ^OB + 



1^R + 

2 dx 



B_ 

2k 



n=l k=l 


-2(-l) n 1 
nB 2 


■^nk^nB 



(-l) n 


'-^nk 


(nB) 2 

-iink^S^ 


(44) 
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that is, the i|r part of the lifting line is one-half of the 

far -wake value plus some effects from the hyperbolic part of the com- 
pressible field. The differentiation in the a direction is made with 
regard for the Fourier development of the saw-tooth function which can 
be recognized in equation (32) at x = 0. For x = 0 it is seen that 
the i in the first brace gives a term 




R.P. 


B_ 

2it ^ 


1 =00 
k=kQ 



-2(-l) n 1 J 

nB 2 


n=l J=Jo 
k=l 


nj 


nk 


JnB 


< 


J nB' 


y nBj 

Pip 

' y nBk 
1 PT 


► (i) 


inBa 


but from equation (17b) this term is just 


1 Bf( p) 

2 2 * 


and, therefore , ' the differentiation of this term cancels one-half of the 
differentiation of the first term in equation (32). In other words, the 
derivative of the above Fourier series with respect to a is to be 
interpreted as the derivative of the sloping part of a saw-tooth func- 
tion of a which the series actually represents, the differentiation 
being performed off the teeth of the saw-tooth function. This choice 
of interpretation is justified by the fact that by so doing the theory 
gives the correct results when degenerated to incompressible flow. 

Since R .P . (^e^ n ^ a ^) contributes nothing at a = g, formal differen- 
tiation of the rest of equation (32) gives 



Therefore , the part of $ f is also one-half of the far-wake value plus some effects from 

the hyperbolic field. 

Next equation (39) is differentiated. The differentiation of the real part along x, for 
x = 0 and a = g, gives 


^(<Pr + <Pb) 


dB -1 


J=°° 

k=k 0 


-2(-l) n 1 


B_ 

^ n=l 116 2 

k=l 




i 


< 

► < 

^nk 

c 


J nB 


¥ nB j 
. P T 

ynBk 
v Pt 


nB 


r (-D n 




nB + 


-i l 


nk 


(nB) 2 

■ iz nkP 2 Ps 2 




(46) 

However, from equation (38), the nB terms in the last brace cancel one-half of the cp^ term 
by the same considerations as for equation (45) • Therefore, equation (46) becomes 



1 x B V" -2(-l) n 1 v T y^Bk 
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and, similarly, 


^-(cpR + 9 b) 


x=0 


a : 


B 


= l 
' 2 8a 


B_ 

2 * 


00 


n=l 



-2(-l) n 

nB 


1 

2 


^nk^nB 



(-l) n (nB) 2 

- 1 iW 2 


(48) 


Equations (44), (45), ( 47 ), and (48) give 0 X and 0 a at the 

lifting line, but they must he expected to include some induction from 
the lifting line itself. This part of the induction must he discarded 
in order to have the Prandtl induced velocities at the lifting line. 

It should be possible to recognize the induction of the lifting line 
from an inspection of equations (44), (45), ( 47 ), and (48). In equa- 
tions (44) and (45), the double summation has coefficients r n ^ which 
depend directly on the circulation distribution T(p) along the lifting 
line; whereas h n k in equations ( 47 ) and (48) depend on h n (p) which 
expresses the field produced by r(p). Therefore, the hyperbolic terms 
in equations (44) and ( 45 ) are now discarded as being caused by the 
lifting line itself. 


Some further explanation of this step is perhaps desirable. It is 
this operation which makes the lifting-line concept applicable 
although the flow has partly supersonic character. The lifting line in 
supersonic flow has infinite wave drag which shows up in the fact that 
the velocities associated with the part of the Reissner potential 
are expressed by series which do not converge (eqs. (44) and (45)). 
Furthermore, if the theory is allowed to degenerate to the case of a 
large number of blades at a large distance from the axis so that the 
potential may be considered independent of radius, the representation 
of an infinite staggered cascade in supersonic flow is obtained, the 
properties of which are already known. Then the hyperbolic terms in 
equations (44) and ( 45 ) are seen to represent just the self-induced 
field of each cascade element by itself. On the other hand, it is the 
self-induced field which must be taken out in order to obtain the 
Prandtl induced velocities. Therefore the indicated rejection is in 
accordance with the lifting-line concept. 
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Substitution of equations (MO and ( 47 ) into equation (12) gives 


*0 


u a = u a . _ A co 

c a i 2n V k=1 


-2(-l) n 


1 

nB 2 




"i^nk - 


(nB) 4 


- U nk3 2 P S 2 


(49) 


and substitution of equations ( 45 ) and (k8) into equation (i+3) gives 

*0 


u t c = u ti + 


_B_ 1 to 
2rt p V 


n=l k=l 


' 2( ' 1)n iw nB ( i ^p')M ^-^4 (50) 


nB 


-1 


■ il nk3 £ 


in which the bars denote the Prandtl induced velocity, the subscripts c 
and i refer to compressible and incompressible, respectively, and the 
terms in P nk have been discarded. 


DISCUSSION 


The Bessel functions in equations (49) and ( 50 ) are those shown in 
figure 2. It is shown subsequently that these functions are distributed 
over the radius in such a way that the first inflection point of the 
Bessel function in the kQth term falls on or very near the sonic radius . 
For the lower k terms the first inflection point is progressively more 
outboard of the sonic cylinder so that at k = 1 it is near the tunnel 
wall with the first extreme (a maximum) of the Bessel function right at 
the tunnel wall. For k = 2, the second extreme (a minimum) is at the 
tunnel wall, and so forth for the succeeding values of k as is seen 
from inspection of the arguments of the Bessel functions (see fig. 3). 
Finally, at k = kQ, which is the hyperbolic term with the largest k, 
the maximum number of oscillations of the Bessel function occur and 
these oscillations are all confined between the sonic cylinder and the 
tunnel wall, because as was pointed out above, the first inflection 
point of the kgth Bessel function falls at or near the sonic cylinder. 

The process of distributing the Bessel functions over the radius 
influences the values of l nk obtained from equation (13). This fact 
brings up the significant point that ^nkQ can approach zero and thus 

produce an effect similar to resonance in a vibrating system; each 
narrow range of conditions for resonance is an extra degree reduced by 
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the square-root operation in equation (13) • An identically zero value 
of Z n kg means that some additional considerations are required. How- 
ever, when later the case of a propeller in free air is taken up, the 
summations over j and k pass to radial integrations in which the 
difficulty of zero Z n appears as a square-root singularity which 

causes no trouble and provides the first hint that a near-zero value 
for Z n k is to be regarded as a necessary feature of the theory due 
to a resonance property of the flow in a tunnel. 

Generation of the coefficients h n k-- The coefficients h n k are 

those in a Fourier-Bessel expansion of the Reissner far-wake radial 
function h n (p). The functions h n (p) have the appearance shown in 
figure 3 . The function h n (p) has a discontinuity in its first deriva- 
tive because, when h n (p) is added to T, the resulting function must 
be smooth in order to make equation (9) give a continuous flow at 
p = p a . Also sketched in figure 3 are the Bessel functions for given 
values of nB used in the expansion of h n (p). These functions have 
their extreme values on the tunnel wall pip as is required by the 

arguments ynB - p. Only the Bessel functions for the hyperbolic terms 
P T 

in the velocity potential are shown, the ones appearing in equations ( 49 ) 
and (50). The value k = 1 corresponds to the Bessel function with the 

y nBk 

first extreme value at p = p T because, from the argument P, 

when p ranges from p = 0 to p T it generates the Bessel function 
up to the first extreme value. The l^th function is the one with first 
inflection point and pg nearly coinciding. Thus, for the case repre- 
sented, k takes on the values k = 1 to k = ko = 4 which correspond 
to hyperbolic terms in the velocity potential. The terms with higher 
numbers of extremes are not shown because they have elliptic flow 
character, and they would merely produce more and more oscillatory por- 
tions inboard of the sonic cylinder p = pg- 

It has been stated that the first inflection point of the kQth 
Bessel function falls at or near the sonic radius. A more precise 
statement, however, would be that the kgth Bessel function has its 
argument equal to its order at or near the sonic radius. However, the 
inflection point has more meaning for visualization purposes . Further- 
more, it can be shown that the argument at the first inflection point 
is nearly equal to the order, if the order is not too small, with incon- 
siderable error as may be seen from figure 2 . This coincidence that 
the kgth Bessel function tends to have an argument equal to the order 
at or near the sonic radius will now be demonstrated and it is to be 
inferred that this coincidence also applies to the first inflection 
point. It will appear that the coincidence of the argument with the 
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order is exact if Z n ^ = 0. For a zero value of Z n kQ, equation (12) 
shows that 


n 02 = v' 
p T ^ ynBkQ 


or, 


^ p = „b (51) 

P T 

which states that the argument of the kgth Bessel function is equal to 
the order nB at the sonic cylinder p = pg if I s exactly zero. 

If ^nkQ no ^ exactly zero, then equation (51) is still nearly 
true. Thus, from equation (12) 

p S = /(nB) 2 + P S 2p2z nk 2 (52) 


which shows that since Zj^ is imaginary, corresponding to hyperbolic 
terms, the argument of the Bessel function at the sonic cylinder is less 
than nB. The reverse would be true if the terms were elliptic so that 
Z n j would be real from equation (12). It appears then that, in fig- 
ure 3, the kQth Bessel function has its argument equal to nB just 
outboard of or at pg, depending on whether Z n k q is merely close to 

zero or identically zero. Further, the hyperbolic terms with k lower 
than kQ have their oscillatory portions more and more outboard of the 
sonic cylinder p = Pg as k decreases from k^ . 

The sonic cylinder and tunnel-wall radius in relation to the first 
appearance of hyperbolic solutions in the velocity po tential.- The 
hyperbolic solutions first enter the velocity potential when the sonic 
cylinder comes inside the tunnel wall. If equation (13) is written in 
the form 
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then since the hyperbolic solutions have imaginary l n ^ the proof of 
this statement is seen to consist in showing that the smallest possible 

I 

value for ^ is unity, which is shown to occur when krj = 1 and 
nB ' 

nB — >00. This proof will place a lower limit on p T for the occurrence 
of imaginary l n \<L’ that is, p T must be greater than p^ for imaginary 
Z nk to occur. Therefore, the solutions cannot change type unless the 
tunnel wall is outside the sonic cylinder. It is obvious that the 

Y T31 

smallest value of must occur for k = 1 because with this value 

nB 

of k the argument of the first extreme must be less than all succeeding 
arguments. It, therefore, remains to show that the smallest value for 

y 1 

- n P — is unity, which follows from a formula on page 14-3 of reference 5 
nB 

that holds when nB is large. 


^nBl _ + O.808618 Yi B ^ 

nB nB 


(55) 


Hence 


nB- 


lim = 1 

„ nB 


(54) 


which completes the proof. 

Diagrammatic point of view of zero l n ^ relation to the 

argument of the Bessel function .- Figure 4 gives a view looking down 
the vertical axis in figure 2 showing the traces of the extremes of 
JnB^y) on the nBjy plane. The slopes of these traces approach unity 
as nB and y become large, a property of the Bessel function. For 
a given propeller operating condition, p T and pg may be calculated; 

Pm 

by plotting in figure 4 the function — nB, it is possible to tell by 

P S 

inspection of the figure how many hyperbolic solutions are present in 
each k-wise summation, and the possibility of a zero value for ^nkQ 

may be appreciated. From equation (51), the condition for zero Z nk ^ 
is 
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(55) 


pij. 


intersects one of the traces of the extremes on an 

Z n kQ will he zero. For 
example , if the dashed line through the origin in figure 4 represents 


If the line — nB 
p S 

integral value of nB, then the corresponding 


pip 


nB for a given operating condition, then Z 


nkQ 


for nB = 8 and 


k = kg = 2 would he zero . 


The dashed line in figure 4 may also he taken as an indication of 
what parts of the Bessel functions JnB(y) are used for the elliptic 
and what parts for the hyperbolic solutions at any given operating con- 
dition. Equation ( 13 ) shows this separation because imaginary Z n k are 

given by equation ( 13 ) so that the solutions are hyperbolic when 



or 



> y nBk 


which means that for hyperbolic solutions, the traces of the extremes 

1 p T 

y nBk raas 't be below the dashed line — nB. Since the slopes of traces 

of the extremes approach unity as nB — the first appearance of 

pip 

hyperbolic solutions occurs when — = 1, because for this case the 
P T p S 

line — nB also has slope unity (shown as the solid straight line 

p S , , 

through the origin in fig. 4). 


The suppression of resonant terms .- The discussion of the resonant 
terms is based on the assumption that in the series expression for the 
induced velocities, or the induced angle of attack at the blade, only 
one zero of Z^v can occur within a finite number of terms in the 


Znk 0 
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series. In the author's opinion, this is true but a rigorous proof is 
lacking. From a practical point of view, the series does not need to 
he computed to a great number of terms, because the higher terms express 
their contributions to the induced angle through oscillatory functions 
of higher and higher frequency so that they eventually become unimportant 
in a practical solution. Thus the question of how many resonant terms 
must be dealt with can be answered for practical purposes by the fact 
that in a series of a reasonable number of terms only one zero of ^nk^ 

can occur at a time so that there is just one resonant term. 


Consider a propeller operating at conditions which maize a particular 
^nkQ exactly zero. With this zero of ^nkQ is associated a resonant 

term or mode in the series expression for the induced angle of attack. 

Now suppose the propeller to be twisted so that the angle of attack 
with respect to the zero-lift angles of the airfoil sections is equal 
to the resonant mode multiplied by some finite constant factor. This 
propeller, according to the theory, is subject to degenerate compressi- 
bility effects , which can be seen without calculation as follows : 


(l) Since infinite induced velocities are physically impossible, 
it is seen from equations (14-9) and (50) that h nk ^ must be zero; that 

is, the lift distribution of the propeller must be such that hnk is 

hnkn n -\) 

zero so that the ratio - — ^ = — appearing in those equations has a 

L nkQ 0 

finite value. 


(2) The only value for hn^/Zni^ which produces zero h n ^ is 

that which makes the resonant induced angle exactly cancel the blade 
twist so that the lift is zero all along the blade. 

Therefore the only solution compatible to both the compressible and 
incompressible components of the induction is that the lift remains 
zero with the indeterminate ratio h n ^/Z nko taking whatever value 

is necessary to make the resonant mode exactly compensate the prescribed 
angle of attack. It seems reasonable to expect a similar suppression of 
the resonant mode when the angle-of-attack distribution is arbitrary. 
Therefore, in compressible flow, the arbitrary prescription of lift 
must always be made with the restriction that h n ] <;Q for the resonant 

mode be zero. The corresponding propeller geometry is seen to be 
indeterminate in that the geometric angle of attack can be altered 
in proportion to the induced angle of the resonant mode without changing 
the prescribed lift distribution. 

The other problem, in which the propeller geometry is given instead 
of the lift distribution, is more realistic. Here, the lift distribu- 
tion must again yield no resonant mode (h n ^ = 0) and is determined 
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according to the requirement that the lifting properties of the blade 
and the induced angles due to all the modes including the resonant mode 
with finite h nko/^nkQ mu st produce the lift distribution by which the 

modes are determined. In practice this solution would be obtained by 
an iteration process. 

Assume a propeller of fixed geometry. The lift distribution, or 
otherwise r(p), is to be determined along the blade. If the resonant 
term is assumed to be zero, such a distribution along with the distribu- 
tion of induced velocities and induced angles can be determined. Sup- 
pose, however, that a resonant term occurs in the induced velocities 
and, correspondingly, in the induced angles. The induced angles are 
then infinitely sensitive to small changes in the lift distribution and 
some special considerations are necessary. 


The lift distribution with the resonant term assumed to be zero 
will in general produce some value of the resonant coefficient (b n ^) r . 

To prevent the introduction of infinite induced angles and thus to make 
the lift distribution consistent with the induced-angle distribution, 
the lift distribution must be so modified as to result in zero values 
of h n ^Q so that bn^Q/^nkQ will be finite • 

hnkrv 

Suppose that — = 1, and compute the corresponding change of 

Z nko 

induced angle Axj_(p) and the accompanying change in circulation dis- 
tribution Ar(p) . From Ar(p) compute Ah n ^Q from equations 10(a) 

and 10(b), and an equation similar to ( 19 ). Now because h n (p) « r(p) 
and Ar(p) a Axi(p), take C so that 


C Ah nkQ = -(h nkQ ) r (56) 

and adjust the induced angles and the circulation by adding C Axq(p) 

^nkQ 

and C AT(p), respectively. This is equivalent to taking = C 

Z nkQ 

where h n ^Q = Znkg = 0* The process is inexact, of course, because of 

the effect of the added AT(p) on the induced velocities, and succes- 
sive approximations are therefore required. 

The above iteration process may be performed mentally on the 
hypothetical propeller defined at the beginning of this section in 
which the propeller was twisted with an angle of attack with respect 
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to the zero-lift angle of the airfoil sections equal to the induced angle 
of attack of the resonant mode multiplied by a constant factor. It will 
be seen that the process gives the correct result, which is zero lift 
over the whole blade. The validity of the iteration process is thus con- 
firmed to the extent that the abovementioned hypothetical propeller 
actually does produce no lift. The latter question appears to require 
experimental verification. This requirement is not surprising when it is 
recalled that the Prandtl lifting-line wing theory, in which lifting 
lines were to be replaceable by finite chords of the same circulation, 
could not be rigorously justified and had to await experimental 
verification. 


VELOCITY POTENTIAL IN FREE AIR 


The velocity potential in free air has a form different from that 
where a tunnel is present; the difference is analogous to. the change 
of a Fourier series into an integral when the interval of expansion 
becomes infinitely great. As far as the induced velocities are con- 
cerned, it is believed that they will be not much different from those 
with the tunnel present. Therefore, the free-air potential will not 
be derived in detail. 

From Hankel's integral (see eq. ( 3 ) on p. ^53 of ref. 6) 


F n (R) 




F n (s)JnB(7s)JnB(7R)s ds 


( 57 ) 


the function cp^ in the Reissner potential as p^, — >■<» can be represented 

by 

<Pr = R-P- r 7fn(7)JnB(7P)e inB(a " x) (-i)d 7 ( 58 ) 

n=l riB J 0 

where F n (R) becomes h n (p) by comparison with equation (36) and the 
function f n (7) is given by 

poo 

fn( 7 ) = / h n (x)J nB ( 7 s)ds (59) 

J 0 
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The function h n (p) used here must suit the boundary conditions of the 
far wake, obtained from equations (lO) by letting 

By analogy with equation (ll) and by using a form of equation (57), 
choose 



(60) 

which is a solution of equation (5) if differentiation under the integral 
is permissible. The physical conditions may now be imposed by methods 
similar to those in the case where a tunnel was present. It will be found 

that the functions C n (7) contain the radical I/7 2 - in the 

V p g 2 

denominator, and this radical is the counterpart of Z n j (see eq. (13)). 
Therefore, the summation over j for the tunnel case passes to an inte- 
gration over 7 where the possible infinity caused by Z n j in the 
denominator becomes in free air a square-root singularity in the inte- 
gration, which should cause no difficulty. 


NUMERICAL EXAMPLE 


A two-blade propeller is assumed to operate at an advance Mach 
number M of 0 . 8 l 93 > with p at the tip equal to unity; therefore, 
the advance ratio V/riD is n, and the circulation distribution is 


r*(p) = \ P \/l - P 2 

6 


(61) 


where 


r* 



r 


(62) 
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is a nondimens ional form for r and v/nD is the conventional parameter 
in which nD is tfte angular velocity in revolutions per unit time multi- 
plied by propeller diameter D. The circulation T* is plotted in fig- 
ure 5 , from equation ( 6 l), along with the corresponding induced quan- 
tities u a ^/V, utj_/V, and ai^, where aq^ is the induced angle of 

attack for incompressible flow expressed in radians . No particular 
significance is to be attached to the factor l /6 in equation ( 6 l) or 
to the precise choice M = 0.8193* 

If equations (49) and ( 50 ) are written as 


where Au a and Auq are the changes due to compressibility, then 


u a c = u ai + 


( 63 ) 


u t c = u ti + Au t 


(64) 



“^nkP^Ps 



( 65 ) 



*#* -X* 

If h n i c , based on r instead of T, is defined as 



( 67 ) 
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the compressibility effects in ratio to the incompressible can be 
determined from equations ( 65 ) and ( 66 ) as follows: 


a c 


Au a /V 

u ai /v 


1 1 

^ s ai /v 


00 lc o 

a nk( p ) 


n=l k=l 


(68) 


■tc 


Au t /V 

u ti /V 


V 

hup ut i /v 



yp) 


( 69 ) 


where 


a nk( P ) = ^ h nk J nBl 


p) 

-17 , - 

\ P T / 

± ^nk p p 


(70) 


. fn \ _ 2 .* T / y nBk \ (nB ) 2 

t n k(p) “ rg hni^nBl ~ — P) r 

113 V p T 


(71) 


The coefficients h nk and -ii nk are tabulated in tables I and II 
for the first six terms in the n-wise summation. These coefficients, 
with the very complete tabulation of the Bessel functions in reference 7 , 
are used to calculate a n k(p) and t nk (p) which are tabulated in 
tables III and IV. 


The ratios 


Aua/V 


and 


Au-t/V 

/V 


u t- 


which express the compressibility 


V v 

effects in ratio to the velocities in incompressible flow are plotted 
against p in figure 6 . The compressibility effect is seen to be very 
small, less than 4 percent. 
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Before discussing the induced angle of attack, some explanation of 
the velocity diagrams in figure 7 may he desirable. Figure 7 (a) shows 
the conventional velocity 'diagram of incompressible-flow propeller 
theory. Figure 7 (b) shows the compressibility effects superposed on 
the incompressible. Since only the part of the diagram showing the 
induced velocities is important, an enlargement of this part is shown 
in figure 7 (c). The lines A-A, B-B, and C-C which are the sides of 
the induced angles may all be drawn parallel to the local helical 
direction in a small-disturbance flow. 

If Axp is the change in induced angle of attack accompanying Au a 
and Au^, then from figure 7 


P Au, 


a 


Au t 


Aai = 


A 


+ p £ 


1/1 


+ p £ 


v /1 + p a 


1 + p c 


Uq 


pa c 


V 


+ t r 


u ts 


(72) 


which, by using equations (68) and ( 69 ), can be put in the form 


where 


A 00 ^ 

a ii ai i 5=1 k=l 


( 73 ) 


Qnk(p) = x ^ ~~ p 2 [p^nk( p ) + £ t nR( P 


( 7*0 


Values of a nk (p) are tabulated in table V. The ratio Ao 1 /a li plotted 

in figure 6 is seen to be very small - , this fact indicates that the induced 
efficiency is reduced only 1 or 2 percent by compressibility, to the 
accuracy provided by six terms in the summation over n. 

In figure 8 , partial sums of o^p) are plotted so that the con- 
tribution of succeeding terms can be seen. There is a shift in the 
convergence with the addition of the fifth term. This shift is caused 
by an unusually large h£ k at n = 5 , kQ = 3- 
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It should be noted that , although the induced angle has been calcu- 
lated for a prescribed circulation distribution, nothing is said about 
the shape of the propeller required to give this circulation distribu- 
tion. At some n, there may be an Z n ^ sufficiently close to zero so 
that a strongly oscillatory induced-angle distribution is finally 
indicated and an unusual blade-angle distribution would result. Appar- 
ently, a degree of arbitrariness is lost in giving arbitrary circulation 
distributions if they are to be related to reasonably shaped propellers. 
More realistic is the view that the propeller shape is prescribed and 
then the performance is calculated. This approach would mean in the 
present example that, after the calculation of the induced angle, its 
effect on the initial circulation distribution would be calculated, 
after which the new distribution would be used for a new induced-angle 
calculation. In the process, the resonant mode would be suppressed, 
and a reasonable induced-angle distribution and a circulation distribu- 
tion somewhat different from that of the prescribed propeller in an 
incompressible flow would result. 


CONCLUDING REMARKS 


The velocity potential for a lifting-line propeller with subsonic 
advance velocity but unrestricted tip speed has been derived in cylin- 
drical harmonics . The theory is for a linear equation of motion and must 
be restricted to advance Mach numbers not too close to one . 

From the velocity potential, the induction at the blade has been 
obtained with a view to preserving a stiict analogy to the induced- 
angle-of -attack concept embodied in Prandtl's lifting-line wing theory 
which has been carried over to propellers by Goldstein and Reissner. 

The Prandtl induced angle of attack becomes representable by the 
basic incompressible results plus some compressibility effects expressed 
in infinite series . Among the terms in the series there is always at 
least one which becomes magnified with indefinite greatness in compari- 
son with the others, an effect resembling resonance in vibration theory. 

A numerical example is provided in which a limited number of terms 
in the series expression for the induced angle of attack have been 
calculated for a propeller with a reasonable lift distribution. The 
compressibility effect appears to be negligibly small, but it must be 
noted that there was no resonant term among the terms calculated. It 
is remarked that the arbitrary prescription of the lift distribution 
cannot be made if resonance is to be considered. 

Although the present calculation indicates that the nonresonant 
terms in the velocity potential produce only a small compressibility 
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effect on the induction at the blade, it is not certain that the resonant 
term, or near-resonant terms, would not make a significant contribution. 

In other wo^ds , if an example were chosen for calculation in which the 
propeller geometry were prescribed instead of the lift distribution, and 
if the operating conditions of the propeller were chosen to make a certain 
convenient term resonant, then, after completing a somewhat laborious 
iteration process, a significant compressibility effect might be observed 
as a result of the direct consideration of the relation of resonance to 
a realistic physical problem. It is the author's opinion that the calcu- 
lation could probably be made with sufficient accuracy by considering 
only the induction of the Reissner potential and the resonant term with 
the nonresonant terms neglected. One might be led further to expect 
such a calculation to explain the experimentally observed dips in the 
lift distributions of propellers operating at tip speeds above a Mach 
number of 1.00; unfortunately, the explanation would not be completely 
clear, however, because there is always the uncertain effect of compressi- 
bility on the airfoil properties in the transonic range which is proba- 
bly a considerable cause of the dips in the lift distribution. 

With regard to experimental verifications there is one check which 
would seem to be satisfactorily definite. It consists in testing a 
propeller in a circular wind tunnel at an operating condition producing 
an exact resonance for a particular term or mode in the series expres- 
sion for the induced angle of attack. If the blade-angle distribution 
for the blade is so designed that the airfoils have, with respect to 
zero lift, an angle-of-attack distribution the same as (or differing 
only by a constant factor from) the induced angle of attack of the 
resonant mode, then, according to the theory, the propeller could 
support no load. The propeller would need to be constructed with some 
accuracy in the blade-angle distribution but the other design particu- 
lars would be arbitrary. In particular, the airfoils might just as 
well be symmetric. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 7, 1953 - 
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TABLE I.- COEFFICIENTS 


\ n 
k \ 

1 

2 

5 

4 

5 

6 

1 

0.010442 

0.008512 

0.005774 

0.005776 

0.002412 

0.001485 

2 



-.000500 

.005185 

.005581 

.005089 

5 





- .005847 

- .000850 


TABLE II.- TABULATION OF -ii nk 


n 

k N. 

1 

2 

5 

4 

5 

6 

1 

5.4966 

7.8164 

12.1428 

16.4777 

20 . 8245 

25.1685 

2 



6.1115 

11.4201 

15-9732 

20.4655 

5 





8.2499 

14.1980 
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TABLE III.- ESSENTIAL PART OF THE DHXICED-AXIAL- VELOCITY CALCULATION a n k(p) 


Nv n 

k 

1 

2 

3 

4 

5 

6 

P = 0.3 

1 

2 

3 

-0.001702 

-0 .000066 

- 0.000003 

.000009 

-0.000006 

0.000017 


1 

k=l 

-.001702 

- .000066 

.000006 

-.000006 

.000017 


P = 0.4 

1 

2 

3 

- 0.002949 

-0.000198 

- 0.000014 

.000043 

-0.000001 

-.000049 

-0.000007 

.000222 

-0.000001 

.000004 

£ 

k=l 

-.002949 

-.000198 

.000029 

-.000050 

.000215 

.000003 

P = 0.5 

1 

2 

3 

- 0.004468 

- 0.000457 

-0.000051 

.000134 

-0 .000006 
-.000229 

-0.000001 

-.000051 

.001384 

-0.000010 

.000037 

£ 

k=l 

-.004468 

-.000457 

.000083 

-.000235 

.001332 

.000027 

p = 0.6 

1 

2 

3 

-0.006192 



-O.OOO881 

r 

-0.000138 

.000308 

-0.000023 

-.000735 

r 

- 0.000004 

-.000226 

.005122 

-0.000001 
- .000062 
.000194 

£ 

k=l 

- .006192 

- .000881 

.000170 

-.000758 

.004892 

.000131 

P = 0.7 

1 

2 

3 

-0.008036 

- 0.001499 

-0.000307 

.OOO56O 

-0 .000066 
- .001758 

- 0.000014 
- .000711 
.012522 

-0.000003 

-.000257 

.000630 

£ 

k=l 

-.008036 

-.001499 

.000253 

-.001822 

.011797 

.000370 

p = 0.8 

1 

2 

3 

- 0.009944 

- 0.002311 

-0 .000593 

.OOO838 

-0.000160 

-.003296 

- 0.000044 

-.OOI67O 

.020977 

-0.000012 

-.000756 

.001352 

la 

-.009944 

-.002311 

.000245 

-.003456 

.019263 

.000584 

p = 0.9 

1 

2 

3 

-O.OH81I 

-0.003298 

-0.001017 

.001040 

-0.000330 

-.004963 

- 0.000109 

-.003040 

.023273 

-0.000036 
- .001662 
.001907 

£ 

k=l 

- .011811 

-.003298 

.000023 

-.005293 

.020124 

.000209 
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TABLE IV.- ESSENTIAL PART OF THE IlffiUCED- TANGENTIAL- VELOC ITY CALCULATION t^p) 


\ n 

k 

1 

2 

3 

4 

5 

6 

rov 

O 

II 

a 

L 

2 

3 

0.001643 

0.000084 

0.000004 

-.000005 

0.000004 

-0.000009 

— 


.001643 

.000084 

-.000001 

.000004 

-.000009 


p = 0.4 

1 

2 

3 

0.002846 

0.000252 

0.000021 

-.000026 

0.000002 

.000035 

0.000006 
- .000122 

0.000001 

-.000002 

*0 

k=l 

.002846 

.000252 

- .000005 

.000037 

-.000116 

-.000001 

P = 0.5 

1 

2 

3 

0.004312 

O.OOO58I 

0.000074 

-.000079 

0.000009 

.000167 

0.000001 

.000042 

-.OOO761 

0.000009 

-.000024 

£ 

k=l 

.004312 

.OOO58I 

-.000005 

.000176 

-.000718 

-.000015 

p = 0.6 

1 

2 
3 

0 .005976 

0.001121 

0.000198 

-.000181 

0.000035 

.000535 

0.000006 

.000188 

-.002818 

0.000001 

.000057 

-.000123 


.005976 

.001121 

.000017 

.000570 

-.002624 

-.000065 

p - 0.7 

1 

2 

3 

0.007756 

0.001909 

0.000442 

-.000330 

0.000102 

.001280 

0.000023 

.000591 

-.006889 

0.000005 

.000237 

-.000398 


.OO7756 

.001909 

.000112 

.001382 

-.006275 

-.000156 

p = 0.8 

1 

2 

3 

0.009598 

0.002942 

O.OOO853 

-.000493 

0 .000247 
.002403 

0.000071 

.001389 

-.011541 

0.000020 

.000697 

-.000855 

£ 

k=l 

.009598 

.002942 

.000360 

.002650 

-.010081 

-.00138 

P = 0.9 

1 

2 

3 

0 .011400 

0.004198 

0.001464 

-.000612 

0.000511 

.003618 

O.OOOI77 

.002529 

-.012804 

0.000060 

.001532 

-.001206 

£ 

k=l 

.011400 

.004198 

.000852 

.004129 

-.010098 

.000386 
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TABLE V.- ESSENTIAL PART OF THE INDUCED- ANGLE -OF -ATTACK CALCULATION o^p) 


n 

k 

1 

2 

3 

4 

5 

6 



P = 0.3 



1 

2 

3 

0.004557 

0.000239 

0.000012 

-.000012 

0.000009 

-0 .000021 


£ 

k=l 

.004557 

.000239 

0 

.000009 

-.000021 



.004557 

.004796 

.004796 

.004805 

.004784 

0.004784 



p = 0.4 



1 

2 

3 

0.005116 

0.000474 

0.000041 

-.000041 

0.000004 

.000058 

.000011 

-.OOOI85 

0.000002 

-.000002 

s 

.005116 

.000474 

0 

.000062 

-.000174 

0 

IT=1 Si 

.005116 

.005590 

.005590 

.005652 

.005478 

.005478 



P » O.5 



1 

2 

3 

0 .005112 

0.000747 

0 .000098 
-.000072 

0.000011 

.000176 

0.000002 

.000046 

-.000664 

0.000010 

-.000024 

Si 

.005112 

.000747 

.000026 

.000187 

-.000615 

- .000014 

£t 

n=l k=l 

.005112 

.005859 

.005885 

.006072 

.005456 

.005442 



P = 0.6 



1 

2 

3 

0.004592 

0.000985 

0.000181 

-.000086 

0.000035 

.000331 

0.000006 

.000131 

-.001194 

0.000001 

.000043 

-.OOOO65 

J: 

k=l 

.004592 

.000985 

.000095 

.000364 

-.001057 

-.000021 

±i 

n=I k=l 

.004592 

.005577 

.005672 

.006036 

.004979 

.004958 



p - 0.7 



1 

2 

3 

O.OO3661 

0.001126 

0.000280 

-.000054 

0.000067 

.000403 

0.000015 
.000233 
- .000722 

0.000004 

.000106 

-.OOOO85 

£ 

k=I 

.003661 

.001126 

.000226 

.000470 

-.000474 

.000025 

±i 

n=l k=l 

.003661 

.004787 

.005013 

.005483 

.005009 

.005034 



P = 0.8 



1 

2 

3 

0.002465 

0.001115 

0 .000361 
.000033 

0.000111 

.000224 

0.000033 

.000244 

.001436 

0.000014 

.000162 

.000008 

Si 

.002465 

.001115 

.000394 

.000335 

.001713 

.000184 

tl 

n=l k=l 

.002465 

.003580 

.003974 

.004309 

.006022 

.006206 


P = 0.9 


1 

2 

3 

0.001125 

0.000937 

0 .000393 
.000141 

0.000150 

-.000247 

O.OOOO55 

.000041 

.003712 

0.000019 

•000114 

•000208 

£ 

k=l 

.001125 

.000937 

.000534 

-.000097 

.003808 

■000341 

t't 

n=l k=l 

.001125 

.002062 

.002596 

.002499 

.006307 

•006648 
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Figure 1 .- Orientation of cylindrical coordinate system. 



Figure 2.- The function J^y) . 
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Figure 3-- Schematic representation of the Bessel functions J^Cy) in 
relation to the functions r(p) and h n (p). 
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Figure 4 .- Traces of the extremes of J^Cy) on the hB,y plane. 


-r 

4 =“ 



Figure 5 • - Circulation distribution and corresponding incompressible 


induced quantities. 



B = 2. 
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Figure 6 .- Compressible induction at the blade. 


A6 


NACA TN 2985 




(a) (b) 


Lift 


ABC 



(c) 


Figure 7 «- Induced velocities at the blade. 
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Figure 8.- Partial sums of a^Cp), the essential quantity in the induced- 

angle-of-attack calculation. 


4 =- 
— J 


$865 NlL VDVM 


